Method for Generating Clock Corrections for a Wide- Area or 
Global Differential GPS System 

FIELD OF THE INVENTION 

[0001] The present invention relates generally to positioning and navigation systems 
that use the Global Positioning System (GPS), and particularly to methods of 
computing satellite clock corrections for wide area or global differential GPS (DGPS) 
networks. 

BACKGROUND OF THE INVENTION 

[0002] GPS uses satellites in space to locate objects on earth. With GPS, signals from 
the satellites arrive at a GPS receiver and are used to determine the position of the 
GPS receiver. Currently, two types of GPS measurements corresponding to each 
correlator channel with a locked GPS satellite signal are available for civilian GPS 
receivers. The two types of GPS measurements are pseudorange, and integrated 
carrier phase for two carrier signals, LI and L2, with frequencies of 1.5754 GHz and 
1.2276 GHz, or wavelengths of 0.1903 m and 0.2442 m, respectively. Some receivers 
output a Doppler measurement which is simply the difference in two consecutive 
integrated carrier phase measurements. The pseudorange measurement (or code 
measurement) is a basic GPS observable that all types of GPS receivers can make. It 
utilizes the C/A or P codes modulated onto the carrier signals. The measurement 
records the apparent time taken for the relevant code to travel from the satellite to the 
receiver, i.e. the time the coded signal left the satellite according to the satellite clock 
minus the time it arrives at the receiver according to the receiver clock. 

[0003] The carrier phase measurement is obtained by integrating a reconstructed 
carrier of the coded signal as it arrives at the receiver. Thus, the carrier phase 
measurement is also a measure of a transit time difference as determined by the time 
the signal left the satellite according to the satellite clock and the time it arrives at the 
receiver according to the receiver clock. However, because an initial number of whole 
cycles in transit between the satellite and the receiver when the receiver starts tracking 
the carrier phase of the signal is usually not exactly known, the transit time difference 



l 



CAl: 344558.1 



may be in error by the periods of a few carrier cycles, i.e. there is a whole-cycle 
ambiguity in the carrier phase measurement. Since the carrier frequencies are much 
higher and their pulses are much closer together than those of the C/A or P codes, the 
carrier phase measurement can be much more accurate than the code measurement. In 
addition, the code measurements are affected much more than the carrier phase 
measurements by the interference of reflected signals with the direct signal. This 
interference also causes the code measurements to be less accurate than the carrier 
phase measurements. 

[0004] With the GPS measurements available, the range or distance between a GPS 
receiver and a satellite is calculated by multiplying a signal's travel time by the speed 
of light. These ranges are usually referred to as pseudoranges (false ranges) because 
the receiver clock generally has a significant time error which causes a common bias 
in the measured range. This common bias from receiver clock error is solved for 
along with the position coordinates of the receiver as part of the normal navigation 
computation. Various other factors can lead to errors or noise in the calculated range, 
including ephemeris error, satellite clock timing error, atmospheric effects, receiver 
noise and multipath error. To eliminate or reduce these errors, differential operations 
are typically used in GPS applications to cancel the noise factors in the pseudorange 
and/or carrier phase measurements resulting from these error sources. Differential 
GPS (DGPS) operations typically involve a base reference GPS receiver, a user GPS 
receiver, and a communication mechanism between the user and reference receivers. 
The reference receiver is placed at a known location and the known position is used to 
generate corrections associated with some or all of the above error factors. The 
corrections are supplied to the user receiver and the user receiver then uses the 
corrections to appropriately correct its computed position. The corrections can be in 
the form of corrections to the reference receiver position determined at the reference 
site or in the form of corrections to the specific GPS satellite clock and/or orbit. 
Corrections to the reference receiver position are not as flexible because, for optimum 
accuracy, they require that the same satellites be observed by the user receiver and the 
reference receiver. 

[00051 The fundamental concept of Differential GPS (DGPS) is to take advantage of 
the spatial and temporal correlations of the errors inherent in the GPS measurements. 
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The GPS satellite clock timing error (clock error), which appears as a bias on the 
pseudorange or carrier phase measurement, is perfectly correlated between the 
reference receiver and the user receiver. So, a DGPS system is capable of completely 
removing the clock error, which typically contributes about several meters of error to 
the user position. 

[0006] The atmospheric effects are due to the GPS signal passing through the charged 
particles of the ionosphere and then through the water vapor in the troposphere. The 
effect of the ionosphere on the GPS signal is usually modeled by an ionospheric 
refraction model and errors in the model contribute to errors in the computed range. 
These errors are strongly correlated between the reference and user receivers over 
short distances between the two receivers, but the correlation diminishes over large 
distances. 

[0007] Refraction of the GPS signal in the troposphere can generally be modeled to 
remove 90 to 95 percent of the tropospheric effects. The residual tropospheric errors 
can be reduced by the use of a DGPS when the user does not roam across large 
distances, because the correlation of tropospheric refraction error typically disappears 
as soon as the user is a few tens of kilometers away from the reference receiver. So 
the use of the DGPS by itself does not produce meaningful reductions of the residual 
tropospheric refraction error. Sometimes, the uncorrected error at the reference 
receiver can even worsen the situation by introducing additional error into the 
computed user position. 

[0008] The ephemeris or satellite orbital error can be modeled as having along-track, 
cross track, and radial error vectors. The satellite orbital error can be reduced by using 
the DGPS system. However, the reduction is somewhat limited because the 
correlation in the satellite orbital error gradually reduces as the separation between the 
reference and user receivers increases. The correlation is largely diminished over 
continental distances. 

[0009] Receiver noise and multipath (reflected signal) effects are generally 
uncorrected between the reference and user receivers. These error effects are 
sometimes amplified by the use of DGPS systems. 
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[0010] To overcome the inaccuracy of the DGPS system in wide-area applications, 
various wide area DGPS (WADGPS) techniques have been developed. The 
WADGPS includes a network of multiple reference stations in communication with a 
computational center or hub. Error corrections are computed at the hub based upon 
the known locations of the reference stations and the measurements taken by them. 
The computed error corrections are then transmitted to users via a communication link 
such as satellite, phone, or radio. 

[0011] In some cases, raw data such as the measurements and positions of the 
reference receivers are supplied to the user receiver(s) rather than the corrections. The 
user receiver can select the data from a particular reference station or form corrections 
from a weighted combination of the data from the multiple reference stations. 

[0012] By using multiple reference stations, WADGPS provides more accurate 
estimates of the error corrections. However, the use of multiple reference stations also 
makes computation of the error corrections more complicated and different error 
factors can alias into each other, destroying the correlations inherent in the GPS 
measurements. 

SUMMARY OF THE INVENTION 

[0013] The present invention includes a method for generating satellite clock 
corrections for a wide area differential GPS (WADGPS) network. In one embodiment 
of the present invention, the WADGPS network includes a plurality of reference 
stations each having a dual-frequency GPS receiver, which facilitates obtaining GPS 
measurements, including pseudorange code measurements and carrier phase 
measurements, on both the LI and L2 carrier signals. With the dual frequency 
measurements, ionospheric refraction effects can be completely removed from the 
GPS measurements taken at the reference stations. Or, alternatively, a linear 
combination of carrier-phase measurements can be formed to match the ionospheric 
refraction effects of the corresponding code measurements. The removal of the 
ionospheric refraction effects or the difference between the ionospheric refraction 
effects on the code measurements and on the corresponding carrier phase 
measurements allows unlimited smoothing of the code measurements with the 
corresponding carrier-phase measurements. This in turn allows a virtual elimination 
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of multipath noise after smoothing for a certain period of time, such as ten minutes or 
more. Satellite orbital errors, because they change much more slowly, can either be 
removed in a separate computational process when the WADGPS network is a global 
DGPS network, or can simply be ignored when the WADGPS network is for a 
continental sized region. The tropospheric refraction effect can be largely removed by 
modeling, and if desired, can be improved by the use of small stochastic adjustments 
included in the computation of the clock correction. After removing the above error 
factors, satellite clock corrections are computed for the individual reference stations. 
An average clock correction is formed thereafter for each of a plurality of satellites by 
taking an average or weighted average of the satellite clock corrections over reference 
stations to which the satellite is visible. 

[0014] Because the method of the present invention computes satellite clock 
corrections after the removal of all other substantial error components, the method 
can be strikingly simple and very robust as compared to conventional methods, which 
often employ Kalman filters to solve for many different parameters simultaneously. 
The disadvantage of these conventional methods is that various error sources can alias 
into one another and destroy the perfect correlation inherent in the satellite clock 
errors. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] Figure 1 is a block diagram of a WADGPS network according to one 
embodiment of the present invention. 

[0016] Figure 2 A is a block diagram of a computer system that serves as an example 
of a main computer system in the WADGPS network according to one embodiment of 
the present invention. 

[0017] Figure 2B is a block diagram of a computer system at a reference station in the 
WADGPS network according to one embodiment of the present invention. 

[0018] Figure 3 A is a flowchart illustrating a method for forming a smoothed 
refraction-corrected code measurement according to one embodiment of the present 
invention. 
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[0019] Figure 3B is a flowchart illustrating a method for smoothing the refraction- 
corrected code measurement according to one embodiment of the present invention. 

[0020] Figure 3C is a flowchart illustrating a method for smoothing the refraction- 
corrected code measurement according to an alternative embodiment of the present 
invention. 

[0021] Figure 4 A is a flowchart illustrating a method for computing satellite clock 
corrections for the WADGPS network according to one embodiment of the present 
invention. 

[0022] Figure 4B is a flowchart illustrating a method for computing satellite clock 
corrections for the WADGPS network according to an alternative embodiment of the 
present invention. 

[0023] Figure 4C is a flowchart illustrating a method for computing satellite clock 
corrections for the WADGPS network according to another alternative embodiment of 
the present invention. 

[0024] Figure 4D is a flowchart illustrating a method for computing satellite clock 
corrections for the WADGPS network according to yet another alternative 
embodiment of the present invention. 

DESCRIPTION OF THE INVENTION 

[0025] FIG. 1 illustrates a wide-area or global DGPS (WADGPS) system 100 
according to one embodiment of the present invention. As shown in FIG. 1, the 
WADGPS network 100 includes a plurality of satellites 1 10 5 a plurality of reference 
stations 120 each having a GPS receiver 122, and communication links (not shown) 
among the satellites 1 10 and the reference stations 120. The reference stations 120 are 
placed at known locations across a wide area 101, such as a continent, for a wide-area 
DGPS system, or across the globe for a global DGPS network. The WADGPS 
network 100 may be utilized by one or more users 140 each having a user GPS 
receiver 142 for positioning and/or navigation purposes. 

[0026] The WADGPS system 100 further includes a main computer system (not 
shown in FIG. 1) that is shared among the reference stations. The main computer 
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system can be located at one of the reference stations or at a different place in or near 
the area 101. FIG. 2 A shows a block diagram of an exemplary computer system 200 
that can be used to as the main computer system. Referring to FIG. 2A, the computer 
system 200 can be a microprocessor-based computer coupled via the internet and/or 
various other communication links to the reference stations 120. Computer system 
200 includes a central processing unit (CPU) 202, memory 210, a multitude of input 
ports 204 and an output port 206, and (optionally) a user interface 208, coupled to 
each other by one or more communication buses 209. The memory 210 may include 
high-speed random access memory and may include nonvolatile mass storage, such as 
one or more magnetic disk storage devices. Memory 210 may include mass storage 
that is remotely located from the central processing unit 202. The memory 210 
preferably stores an operating system 212, a database 214, and GPS application 
procedures 216, including procedures 218 for implementing methods of generating 
satellite clock corrections in the embodiments of the present invention, as described in 
more detail below. The operating system 212 and application programs and 
procedures 216 and 218 stored in memory 210 are for execution by the CPU 202 of 
the computer system 200. The memory 210 preferably also stores data structures used 
during execution of the GPS application procedures 216 and 218, including the 
smoothed refraction corrected code measurements and mean receiver clock errors 
associated with individual reference stations as described below, as well as other data 
structures discussed in this document. 

[0027] The operating system 212 may be, but is not limited to, an embedded 
operating system, UNIX, Solaris, or Windows 95, 98, NT 4.0, 2000 or XP. More 
generally, operating system 212 has procedures and instructions for communicating, 
processing, accessing and storing data, and for performing other basic operations. 

[0028] The input ports 204 are for receiving data from the reference stations 120, and 
output port 206 is used for outputting data and/or calculation results. Data and 
calculation results may also be shown on a display device of the user interface 208. In 
one embodiments of the present invention, GPS measurement data taken at the 
reference stations are sent to the input ports 204 of the computer system 200 for 
processing by the computer system 200 according the procedures 218, as described in 
more detail below. Calculation results are output through the output port 206 of the 
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computer system 200 and are sent via redundant communication links to an upload 
station among the reference stations 220. The upload station sends the correction data 
to one or more of the satellites 210 for broadcast to the user 240. 

[0029] The computer system 200 usually includes a very high-speed processor since 
it may need to handle raw GPS measurement data from a multitude of GPS receivers. 
A continental WADGPS system usually has about 3 to 10 reference receivers and a 
global WADGPS system usually has about 20 to 100 reference receivers feeding data 
to the computer system 200. 

[0030] In some applications, each of the reference stations 120 may include a 
computer system 124 coupled to the GPS receiver 122. As shown in FIG. 2B, a 
computer system 124 coupled to the GPS receiver 122 at a reference station 120 
includes a central processing unit (CPU) 126, memory 128, an input port 134 and an 
output port 136, and (optionally) a user interface 138, coupled to each other by one or 
more communication buses 129. The memory 128 may include high-speed random 
access memory and may include nonvolatile mass storage, such as one or more 
magnetic disk storage devices. The memory 128 preferably stores an operating 
system 131, a database 133, and GPS application procedures 135. The GPS 
application procedures may include procedures 137 for implementing methods of 
forming smoothed refraction corrected code measurements in the embodiments of the 
present invention, as described in more detail below. The operating system 131 and 
application programs and procedures 135 and 137 stored in memory 128 are for 
execution by the CPU 126 of the computer system 124. The memory 128 preferably 
also stores data structures used during execution of the GPS application procedures 
135 and 137, including GPS pseudorange and carrier-phase measurements 139, as 
well as other data structures discussed in this document. 

[0031] The input port 134 is for receiving data from the GPS receiver 122, and the 
output port 136 is used for outputting data and/or calculation results to the main 
computer system 200 of the WADGPS system. Data and calculation results may also 
be shown on a display device of the user interface 138. In one embodiment of the 
present invention, the CPU 126, the memory 128 and the input port 134 of the 
computer system 124 are integrated with the GPS receiver 122 into a single device, 
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within a single housing, as shown by the dashed line in FIG. 2B. However, such 
integration is not required to carry out the methods of the present invention. 

[0032] In one embodiment of the present invention, satellite clock corrections are 
computed using smoothed refraction-corrected code measurements associated with 
one or more reference stations. For each of the one or more reference stations, the 
smoothed refraction-corrected code measurements are formed at each measurement 
epoch from the pseudorange and carrier phase measurements taken by the GPS 
receiver at the reference station. Many GPS receivers make both a C/A-code 
measurement and a P-code measurement on the LI or L2 frequency, and either of the 
C/A or P-code measurements can be used as the LI or L2 code measurement. 
However, since small biases exist between the two measurements, whichever of the 
two is used in the reference receiver, the same should also be used for the equivalent 
process in the user GPS receiver(s). In the following discussion, for each satellite 
visible at the reference station, the LI and L2 frequencies will be designated as fi and 
f 2 , respectively, the raw pseudorange code measurements on the LI and L2 
frequencies at a measurement epoch will be designated as Pi and P 2 , respectively, and 
the raw carrier phase measurement on the LI and L2 frequencies at the same 
measurement epoch will be designated as <f>\ and <fc. 

[0033] FIG. 3A illustrates a method 300 for obtaining a smoothed, refraction- 
corrected code measurement from GPS data taken at a reference station according to 
one embodiment of the present invention. As shown in FIG. 3 A, method 300 includes 
step 310 in which the raw code and carrier-phase measurements for each satellite 
visible at the reference station are first refraction corrected to obtain refraction 
corrected code and carrier-phase measurements. The refraction-corrected code 
measurement, designated as Prc, is formed as follows: 

P»c = .J'^ P* ~ ,/ 2 \, P 2 =P* -1.5457(P, -P 2 ) 



(/,'-//) (/,"//) . (1) 



The refraction-corrected carrier-phase measurement, designated as L RC , is formed in 
similarly as follows: 
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where Li and L 2 are the carrier-phase measurements scaled by the wavelengths of the 
LI and L2 signals, respectively, and each includes an approximate whole-cycle 
ambiguity value that has been added to cause the scaled carrier-phase measurement to 
be close to the same value as the corresponding code measurement. Thus, 



where the whole-cycle values of Ni and N 2 have been initialized at the start of carrier- 
phase tracking to give values that are within one carrier wavelength of the 
corresponding code measurements so as to keep the differences between the scaled 
carrier-phase measurements and the corresponding code measurements small. From 
the form of equation (2), it is noted that the refraction corrected carrier-phase 
measurement includes a whole-cycle ambiguity with a wavelength determined by the 
sum of the fi and f 2 frequencies, which is approximately 0.1070 meters. 

[0034] Because the ionospheric effects have been removed from both the code and 
carrier-phase measurements by the operations performed in step 310 according to 
equations (1) - (4) and the effects of satellite clock and orbit errors on the 
pseudorange and carrier-phase measurements are the same, the values of Prc and Lrc 
obtained in step 310 should be almost identical except for the possible whole-cycle 
ambiguity associated with the carrier-phase measurements and the higher multipath 
noise in the code measurements. This allows the formation of a carrier-phase 
smoothed code measurement, which approaches the small measurement noise of the 
carrier-phase but without the associated whole-cycle ambiguity. 

[0035] Thus, method 300 further includes step 320 in which the refraction corrected 
code measurement is smoothed by the corresponding refraction corrected carrier- 
phase measurement. In one embodiment of the present invention, as shown in FIG. 
3B, step 320 is performed in two substeps, substep 322 and substep 324. In substep 
322, for each satellite visible at the reference station, an offset between the refraction 
corrected code measurement and the refraction corrected carrier-phase measurement 



(3) 



L 2 =(<p 2 +N 2 )Z 



(4) 
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at each of a series of measurement epochs is computed and an expanding average is 
taken to form a smoothed offset Oi as follows: 



O|=0H+(^C-Ac-0 H )/tf, 



(5) 



where i is used to designate the current measurement epoch, and the value of t\ is 
equal to i until a maximum value of averaging is attained. For example, if the carrier- 
phase measurement is assumed to have only 1/1 00 th of the noise of the code 
measurement, the value of "if would be limited to 100 squared or 10,000. 

[0036] In substep 324, the smoothed, refraction-corrected code measurement, S, is 
obtained by adding the refraction corrected carrier-phase measurement for the current 
measurement epoch to the smoothed offset, so that 



[0037] In an alternative embodiment, the two-step process for performing step 320 as 
described above can be combined into a single-step, so that: 



This single step process projects the measurement ahead using the change in the 
refraction corrected carrier-phase measurement and then averages the difference 
between that projection and the code measurement. However, the two-step process 
has an advantage in that the constancy of the offset value can be monitored and, if 
desired, a threshold change in the value can be used to edit out erroneous 
measurements. 

[0038] FIG. 3C illustrates another method 350 for obtaining the smoothed refraction- 
corrected code measurements from GPS data obtained at a reference station, 
according to an alternative embodiment of the present invention. As shown in FIG. 
3C, method 350 includes step 360 in which a linear combination of Li and L 2 for each 



S,-0,+L, 



(6) 




(7) 



where 



7-1 
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satellite is formed to match the ionospheric refraction effects on each of the code 
measurements Pi and P 2 . The carrier-phase combination that matches the ionospheric 
refraction effect on the Pi code measurement is designated as M\ and is formed as 
follows: 

fx ~fi fx -fi . (8) 

The carrier-phase combination that matches the ionospheric refraction effect on the P 2 
code measurement is designated as M 2 and is formed as follows: 

M 2 =^jtL x -^^L 2 z 5.09Z, -4.09Z 2 

fx ~fz fx ~fi . (9) 

[0039] Method 350 further includes step 370 in which, smoothed code measurements 
can be formed in a fashion parallel to equations (5) and (6) as follows: 

O,, =O f _ u + {P U -M.j-O^/rj, (10) 

or, alternatively, parallel to equation (7) as follows: 

^=£,,+-(^,,-5,,), (12) 

where: 

where subscript j is used to designate the measurements at the two different 
frequencies so that j = 1 or 2. 

[0040] Method 350 further includes step 380 in which the final smoothed refraction- 
corrected code measurement is obtained via the following equation: 

{fx - ft) (fx ~fl) (13) 
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[0041] The value of the smoothed refraction-corrected code measurement obtained 
using equation (7) from the first method is mathematically equivalent to the value of 
the smoothed refraction-corrected code measurement obtained using equation (13) 
from the second method. However, since the refraction correction process amplifies 
noise, any edit threshold used for the change in the refraction corrected offset in 
equation (5) must be larger than a corresponding edit threshold of equation (10). 

[0042] To compute satellite clock corrections, the smoothed refraction corrected code 
measurement is used to compute a residual for the corresponding GPS measurements 
as follows: 

R[ = T k ! -S l k (14) 

where R l k represents the residual, subscript k and superscript / represent the particular 
reference station and the particular satellite, respectively, for which the residual R l k is 
formed, and T l k represents a theoretical range between the reference station k and the 
satellite /. Subscript i representing the current measurement epoch is suppressed in 
Equation (14) and in all of the following discussions. The theoretical range T l k may be 
calculated conventionally using the known position of the reference station. The 
calculation of the theoretical range may include an adjustment or corrections for the 
tropospheric refraction effects computed based on a theoretical model for the 
tropospheric refraction effects. The theoretical range may also include adjustment or 
corrections for the satellite orbital errors. These corrections may be computed in a 
separate or off-line module. They may also be computed less frequently than the 
computation of clock corrections. 

[0043] In one embodiment of the present invention, to generate the satellite clock 
corrections, a master station is selected among a group of reference stations including 
some or all of the reference stations in the WADGPS network. The master reference 
station can be selected based on a judgment that its clock is the most accurate among 
the group of reference stations, or a determination that it has measurements available 
from the largest number of satellites among the group of reference stations, or some 
combined function of these or other factors, or even arbitrarily. The reference stations 



13 



CAl: 344558.1 



other than the master reference station are referred to below as local reference 
stations. 

[0044] FIG. 4A illustrates a method 400 for generating satellite clock corrections for 
the WADGPS network 200 according to one embodiment of the present invention. As 
shown in FIG. 4A, method 400 includes step 410 in which satellite clock corrections 
are computed at the current measurement epoch i for the master reference station. In 
one embodiment of the present invention, for the master reference station, a master 
clock correction for a satellite visible at the master reference station is computed as 
follows: 

rt=R l 0 -M 0 , (15) 

where subscript 0 is used to represent the master reference station, r l Q represents the 
master clock correction for satellite / as computed, and M 0 represents a master mean 
receiver clock error. Note that in the above equation and in some of the following 
equations, the satellite clock corrections are expressed in terms of distance for 
application as direct corrections to the measured ranges. The master mean receiver 
clock error is computed as follows: 

M 0 =— 5> 0 \ (16) 

where the summation is over a group of satellites visible at the master reference 
station, X represents one of the group of satellites, and m 0 represents the number of 
satellites in the group of satellites. 

[0045] Method 400 further includes step 420 in which satellite clock corrections are 
computed at the current measurement epoch i for one or more local reference stations 
each having at least one common satellite with the master reference station. A 
reference station has a common satellite with another reference station, or the satellite 
is common to both of the reference stations, when both of the reference stations have 
GPS measurements available from the satellite to form the smoothed refraction- 
corrected code measurements for the satellite at the current measurement epoch as 
described above. 
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[0046] For a local reference station, a local clock correction for a satellite common to 
the local reference station and the master reference station is computed as follows: 

rl=R' k -M k , (17) 

where r[ represents the local clock correction for local reference station k for satellite 
/, and M k represents a local mean receiver clock error. The local mean receiver clock 
error is computed as follows: 

M k =— -r?) = M 0 + AM A (18) 

m k x 

where the summation is over a group of satellites visible at both the master reference 
station and the local reference station k, X represents one of this group of satellites, m k 
represents the number of satellites in this group of satellites, and 

AM, = -Ri) =— £(** -r 0 *)-M 0 

represents the difference between the local mean clock error for the local reference 
station k and master mean receiver clock error. The summation in Equation (18) is 
over satellites common to the master reference station and the reference station k, and 
m k represents the number of such satellites. Once M k is obtained using satellites 
common to the master reference station and the local reference station in equation 
(18), local clock correction for other satellites visible at the local reference station and 
not to the master reference station can be computed by computing equation (17) for 
those satellites. 

[0047] After obtaining satellite clock corrections for the individual reference stations 
in step 420, method 400 further includes step 430 in which an average clock 
correction is computed for every satellite visible either at the master reference station 
or visible at some or all of the local reference stations having one or more common 
satellites with the master reference station. In one embodiment of the present 
invention, an average clock correction is computed for a satellite in step 430 by 
averaging the clock corrections computed for the satellite in step 420, as expressed by 
the following equation: 
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p k 

where the summation is over reference stations to which satellite / is visible and p is 
the number of such reference stations. 

[0048] FIG. 4B illustrates a method 450 for generating satellite clock corrections for 
the WADGPS network 200 according to an alternative embodiment of the present 
invention. In method 450, the satellite clock corrections are computed iteratively. In 
each iteration, Equations (17) and (18) are computed for each of a subgroup of local 
reference stations, and the average in Equation (19) is over the master reference 
station and the subgroup of local reference stations. The subgroup of local reference 
stations originally includes no reference station and a local reference station is added 
to the subgroup of local reference stations in each iteration. 

[0049] As shown in FIG. 4B, method 450 includes step 455 in which satellite clock 
corrections are computed for the master reference station according to Equations (15) 
and (16). Method 450 then proceeds to step 460 in which a first local reference 
station is added to the subgroup of local reference stations. The first local reference 
station is preferably a local reference station having the most number of common 
satellites with the master reference station. Method 450 then proceeds to step 470 in 
which satellite clock corrections are computed for the first local reference station 
using Equations (15) through (18), and to step 480 in which average clock corrections 
are computed by averaging the satellite clock corrections computed for the master 
reference station and the first local reference station using Equation (19). 

[0050] If there are other local reference stations having common satellites with the 
master reference station, method 450 goes back to step 460 to add one of these other 
local reference stations to the subgroup of local reference stations. In general, a 
reference station sharing the next largest number of common satellites with the master 
reference station is selected. Then, method 450 proceeds to perform the calculations 
in step 470 and 480 for the new subgroup using equations (17) through (19) except 
that the term r 0 ' in Equation (18) is replaced by c l computed using Equation (19) in 
step 470 in the previous iteration. Thus, the iteration continues until all of the local 
reference stations sharing common satellites with the master reference station have 
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been included in the subgroup of local reference stations and average clock 
corrections are calculated considering the subgroup of local reference stations. 

[00511 The iterative computation of Equations (17) through (19) is especially useful 
for a global DGPS network because a problem can occur in a global DGPS network 
when there are few if any common satellites visible at both the master reference site 
and the particular reference site for which satellite clock corrections are being 
computed. This can cause equation (19) to be either poorly defined or even undefined, 
and the problem can be overcome by the iteration process in method 450. 

[0052] Method 400 or 450 for generating satellite clock corrections is simple and 
effective. Various modifications to the method can be made to improve the accuracy 
of the computed satellite clock corrections without departing from the sprit and scope 
of the present invention. For example, method 400 or 450 can be modified so that a 
weighted combination of the satellite clock corrections computed for the individual 
reference stations is used on the right hand side of Equation (19) to form the average 
clock correction. The individual reference stations can be weighted in equation (19) to 
reflect the angular positions of the satellites visible at the reference stations, or an 
apparent noise level of the GPS receiver at each reference station. 

[0053] Method 400 or 450 may also be modified to allow the mean clock error M k to 
change only slowly from one measurement epoch to the next measurement epoch by 
using a conventional smoothing filter. If the mean clock error Mk is allowed to change 
by a large discrete jump, it can cause step changes in the clock corrections which then 
alias into the user receiver clock solution. This can cause a navigation error if the user 
GPS receiver employs a navigation algorithm that does not allow step changes in the 
receiver clock solution. 

[0054] Method 400 or 450 may also be modified by adding a step 425 between steps 
420 and 430 in method 400, as shown in FIG. 4C, or step 475 between steps 470 and 
480 in method 450, as shown in FIG. 4D. In step 425 or 475, the computed values of 
the clock corrections for the individual reference stations are compared with each 
other, and if the clock correction for one reference station is found to differ 
substantially from those for the other reference stations, it can be edited or replaced 
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with a substitute value before the average clock corrections are computed in step 430 
or 470. 

[0055] Method 400 or 450 may also be modified by forming a weighted combination 
of R l 0 at the right hand side of Equation (16) and/or a weighted combination of 
(/^ _ r J) a t the right hand side of Equation (18) according to the elevation angles of 
the satellites over which the summations in the Equations are taken. Also, if the value 
of (R l k - r l Q ) for one satellite in equation (16) differs substantially from those obtained 
for the other satellites, it can be edited or replaced with an appropriate substitute 
value. 

[0056] An evaluation of the spread of values of {R\ - r l Q ) entering into the 
summation on the right hand side of equation (18) can also be used to adjust on a site 
by site basis the model used to compute the tropospheric refraction effect. A small 
stochastic adjustment could improve the response of the computed clock corrections 
to local weather conditions. This can be accomplished by a procedure which adjusts 
an overhead tropospheric refraction coefficient such that difference between 
{R l k - r J ) a t different reference sites is minimized. In such a process the data from 
satellites with a low elevation may be weighted higher since they are more sensitive to 
the tropospheric refraction effects. 

[0057] Methods, 300, 350, 400, and 450 as described above may be carried out by the 
main computer system 200 of the WADGPS system with the reference stations 
supplying the GPS pseudorange and carrier phase measurements. Alternatively, the 
smoothed refraction corrected code measurements and the residuals may be computed 
at each of the individual reference stations by the computer system 124 coupled to the 
GPS receiver at the reference station. The computation results are then transmitted to 
the main computer system 200 for further processing according to the embodiments of 
the present invention. 

[0058] Thus, the present invention provides a simple and effective method for 
generating GPS satellite clock corrections for a WADGPS network. The computed 
satellite clock corrections are not affected by the ionospheric refraction effects 
because they are removed by the use of dual-frequency measurements. The satellite 
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orbital errors, because they change much more slowly, can either be removed in a 
separate computational process when the WADGPS is a global DGPS network, or 
may simply be ignored when the WADGPS network is for a continental sized region. 
The tropospheric refraction effect can be largely removed by modeling and if desired 
can be improved by the use of small stochastic adjustments included with the clock 
correction computation. The multipath effects are removed from the code 
measurements by smoothing with the corresponding carrier-phase measurements. 
After removing these error factors, the method then provides a simple but very 
accurate computation of the satellite clock corrections for either continental-sized or 
global GPS networks. 
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